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We discuss the Pauli Hamiltonian including the spin-orbit interaction within an U(l) X SU(2) gauge theory 
interpretation, where the gauge symmetry appears to be broken. This interpretation offers new insight into 
the problem of spin currents in the condensed matter environment, and can be extended to Rashba and 
Dresselhaus spin-orbit interactions. We present a few outcomes of the present formulation: i) it automatically 
leads to zero spin conductivity, in contrast to predictions of Gauge symmetric treatments, ii) a topological 
quantization condition leading to voltage quantization follows, and iii) spin interferometers can be conceived 
in which, starting from a arbitrary incoming unpolarized spinor, it is always possible to construct a perfect 
spin filtering condition. 
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electron studies, synchrotron-source X-ray scattering, etc.) 

85.75.-d Magnetoelectronics; spintronics: devices exploiting spin polarized transport or integrated magnetic 
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1. Introduction 

A reformulation of the spin-orbit (SO) coupling Hamiltonian in terms of non-Abelian gauge 
fields [Q] was explicitly given in ref. where the SO interaction is presented as a SU{2) x 

U(l) gauge theory. As the Yang- Mills gauge theory is well understood and is the underpinning of 
well established theory, enormous insight can be brought upon new problems. Such gauge point of 
view, in more general terms, has been known for some time [ [HI [5] [H] • This formulation is very 
revealing, since the consistent gauge structure of the theory becomes obvious and the physics of spin 
currents, persistent currents and color diamagnetism [ UOj can be understood in a manner analogous 
to the well known U(l) gauge theories. A consistent SU{2) x (7(1) gauge approach was presented 
in reference [[HE] where it was found that for the Pauli type Hamiltonians (including Rashba and 
two-dimensional reductions of the Dresselhaus Hamiltonian), Gauge Symmetry Breaking (GSB) is 
necessarily built into the theory and leads to vanishing of the spin conductivity in constant electric 
fields [E]- In addition, the Yang Mills interpretation of the Rashba and Dresselhaus SO interactions 
renders the associated gauge fields real, with topological consequences analogous to the Aharonov 
Casher effect [HE]. 

The Rashba and Dresselhaus SO interactions arise in materials which lack either structural or 
bulk inversion symmetry, respectively [ [HI [T2] [13] . These two kinds of interactions have recently 
been given a great deal of attention due to their potential role in the generation and manipulation 
of spin polarized currents, spin filters [I14j. spin accumulation [[15], and spin optics [I16j. 
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2. Spin-orbit interaction in semiconductors 

Spin-orbit interaction may be defined as the interaction of the spin magnetic moment of the 
electron, fi = —g^s = with the magnetic field produced by all external moving charges 

in the electron rest frame. Here, s = ^frcr is the spin of the electron with cr the vector of Pauli 
matrices, g — 2 the Lande factor, m the electron rest mass and — e the electron electric charge. 
Bold face symbols denote matrices in spin space. In the rest frame of the external charges, the SO 
interaction, properly corrected to take into account Thomas precession, is usually described by the 
Pauli term, ~ 2 \ (px VV), where V is the potential energy of the electron in the presence of 
the external charges. For a spherically symmetric potential, W = f(r)r and the SO interaction 
is proportional to s ■ L, hence its name. 

In a semi-conductor, the electron rest mass and the Lande factor must be replaced by their 
effective counterparts, and the spin-orbit interaction can thus be considerably enhanced. Moreover, 
an appropriate basis in a crystalline solid is given in terms of the Bloch wave functions, ip(f) = 
uge ifcr , and the SO interaction must be calculated in this basis. This is generally a hard task 
and one usually uses phenomenological expressions compatible with the crystal symmetries. When 
there is bulk inversion asymmetry (BIA), we have the Dresselhaus expression in 3d systems [I17j. 

H D ,3d = const. k x {k 2 y - k 2 z )a x + c.p. (1) 

where c.p. stands for cyclic permutations. When the electrons are confined in two dimensions, the 
expectation value along the third dimension should be considered, with (k z ) ~ 0, (fcf) ~ (n/a) 2 , a 
being the typical confinment length in the z— direction. The Dresselhaus SO interaction thus takes 
the simple form 

HD,2d = f3(k x cr x - ky(Ty) (2) 

where we neglect cubic terms in k. If the confining potential is not symmetric i.e. in the case of 
space inversion asymmetry (SI A), there is another term which directly follows from the expression 
A ^i c i & ■ (k x VV) with only the (V Z V) contribution. This term is known as the Rashba SO 
interaction, 

H R £d = a(k y cr x - k x (T y ). (3) 

Note that the Rashba SO amplitude can be tuned experimentally using a gate voltage, since the 
coefficient a is proportional to the electric field. For more details on SO interactions in semi- 
conductors, see Refs. [ H3l H7] . 



3. Non-Abelian gauge field theory approach 



Let us first consider the Pauli SO interaction. Rashba and Dresselhaus SO interactions will be 
treated later. Neglecting electron-electron interactions, we start by considering the Pauli Hamilto- 
nian to order v/c, acting on two-component spinors, 



H = 



(p - eAf 



2m 



+ V -e<j> 



»2x2 



P 



2m 



8to 3 c 2 
ehd ■ (p — eA) x E 
Am 2 c 2 



+ 



^h 2 



8m 2 c' 



;V ■ E 



12x2 



(4) 



where the first term in the first line corresponds to the usual Schrodinger equation including the 
kinetic energy with a mininal coupling to the electromagnetic field, the substrate potential denoted 
by V, that can be assumed periodic, and a scalar potential contribution. The second term in the 
first line describes the first relativistic correction to the kinetic energy and the Darwin term, where 
E is the electric field and c the speed of light. These first two terms are proportional to the 2x2 
identity matrix in spin space. The second line comprises explicitly spin-dependent terms, first the 
Zeeman interaction where B is the magnetic field and <x is the Pauli matrix vector and the second 
term is the spin-orbit interaction, now written with the minimal coupling to the gauge vector. We 
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have assumed a static potential so that the rotor of the electric field is absent and the spin-orbit 
interaction is limited to the term mentioned here. In the following, we absorb the spin-independent 
one-body interactions (second term of first line in Eq. [4]) in the potential, and we rewrite the non 
relativistic kinetic energy plus the spin-orbit interaction as 



K.E. 



s.o. = -L 

2m 
1 

~ 2m 



(p-eA) 2 l 2x2 -2{p-eA)l 2 



x2 



(p-eA) 2 l 2 



eh _ ^2 
a x E 



Amc 2 



Amc 2 
e 2 ft 2 



ax E 

a x E\ 2 . 



32to 3 c 4 



(5) 



This suggests an SU(2) x U(l) form described by non-Abelian W^ a and ordinary Abelian gauge 
fields = (A , A 1 ) = (cf>/c,A) with E = -V0 - d t A and B = V x A, and we can rewrite the 
Hamiltonian, following Jin, Li and Zhang [ I18j as 



H = 



l r 

2m 



{p-eA)l 2x2 -gW a r a 



-— g 2 W a W a l 2x2 - gcW 0a T a + {V - ecA°)l 2x2l 

8m 



(6) 



where the Zeeman interaction is written as the time component of the non-Abelian gauge field 



HH] -gcW° 

gW la T a = 



-a • B, while the space components of this SU(2) connection are defined by 



(eh/2mc )ei a jE J T a , or explicitly, 



gW 1 
gW 2 
gW 3 



eft 

2mc 2 
eft. 

2mc 2 
eft. 

2mc 2 



(E z u y - E y u z ), 
(-E z u x + E x u z ), 
(E y u x - E x u y ). 



(7) 
(8) 
(9) 



with u x , u y and u z unit vectors in the x 



y— and z— directions. The 2x2 matrices r a are the 
symmetry generators for 5/7(2) obeying the commutation relation [r a ,r b ] = ie a b c T c , e a bc being 
the totally antisymmetric tensor. The coupling constant g is fixed by the combination gW^ a which 
has the dimensions of 2 2 E, and we choose g = ft,. The relation between the spin operator and the 
corresponding generators is %T a = s° and the spin is s = 1/2. The superscripts of the beginning 
of the Latin alphabet, a, b, c, ... refer to the internal spin degrees of freedom for which we use 
the convention of summation when they are repeated, while Greek indices /x, v, ... correspond to 
space-time components and run from to 3, the time component corresponding to and the space 
components being also denoted as Latin indices from the middle of the alphabet, i, j, fc, . . . 

This formulation differs from that of ref. [ 1181110] most importantly in the term quadratic in W a 
in equation[BJ The purpose of writing the second term here, as a function of the SU(2) connection 
is to evidence gauge symmetry breaking (GSB) in this Hamiltonian. This observation has important 
consequences in the physical interpretation of the resulting Yang-Mills fields and is the reason why 
the Yang-Mills fields themselves are observable quantities, where as in a gauge symmetric theory 
they would be gauge dependent [ HH [20] . 

In order to generate the Noether currents in a canonical fashion, one must formulate the ap- 
propriate Lagrangian for the corresponding equations of motion. The non-relativistic Lagrangian 
density we seek is (now omitting identity matrices) 



8m 



2 

~4m ^ 4m°^ U ' 



(10) 
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where * is a Pauli spinor * = ( ^ T ) , G°„ = d^W? - d v W% - e abc W b W^ and = d^A v - d v A^ 



#-*+W 6 W 6 * is due to gauge 

8m & & 



are the SU{2) and U(X) field tensors respectively. The new term 
symmetry breaking. The covariant derivative is then of the form —ifiD = —itiV — eA — gW a T a . 
This form of the covariant derivative determines the well known U(l) coupling constant e/h and 
the SU(2) coupling constant for this theory is g/h. A gauge transformation of the Lagrangian 
density would leave it unchanged up to a divergence at the condition that the Coulomb gauge for 
the SU{2) connection is satisfied, V • W a = 0. 

The Hamiltonian in Eq.[5]is derived from the corresponding Lagrange equations for the matter 
fields 'P. The equations of motion of the Yang-Mills fields in the presence of currents (generalized 
Maxwell equations) follow from Euler-Lagrange equations with respect to variations of the gauge 
fields, 



dC 



dC 



(11) 



and since the non-Abelian field tensor 



The l.h.s. of this equation gives ^ 3(3 3 ff°) 
Qfiva j g an tisymmetric with respect to space-time indices, it is natural to introduce a conserved 
current J va — -§^s- The conserved current J va [[!§] is the full spin current carried both by matter 

and radiation i.e. J a = ^ atter + ^Radiation- This full current is an observable physical quantity, 
since the gauge is fixed in the present case (in contradistinction with ordinary SU(2) gauge theory 
where it is gauge-dependent). The spatial (spin current) and time (magnetization) components of 
the current density then follow as 



r = f 

2m 



V4m / 

+ -e ab cW b (G» 
m 



+ G vyc u y + G uzc u z ), 



(12) 



and the spin polarization 



J Qa = *V a * + ^-s abc W b G ]0c . 

m J 

Three terms can be distinguished in the spin current, 
i) The first term has the canonical form expected for the material current namely 



(13) 



Matter = {g / {r a V + VT a )* , 



(14) 



where v % = (l/iK)\r l , H]. 



ii) The second term comes from the gauge symmetry breaking contribution. It depends on both 

matter and radiation. 

iii) Finally, the third term is the canonical radiative contribution originating from the derivative 

with respect to the gauge potential of the non-Abelian contribution of the field tensor G a ^ v . 
Note that this last term would not arise in an Abelian theory (e.g. in C/(l), since the photon 
does not carry any electric charge). 

The first and last terms were described in ref. [ [IB] as taken from an apparently gauge symmetric 
form. The magnetization term has both a material contribution (the first term) and a radiative 
contribution as both matter and radiation carry angular momentum. We emphasize that, the extent 
to which gauge symmetry is broken depends on the choice of the electric field. If only the E z ^ 0, 
then one allows gauge transformations that leave W% — —W\ invariant. This is analogous to the 
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(r a <i') t W- (V*) f (x a *) -&(2-W a ^V ( 15 ) 



remnant -Z2 group after U(l) GSB in superconductors and a similar situation in the electro- weak 
GSB mechanism [ITS]. 

Since the gauge theory considered here is non-Abelian, in the most general case the correspond- 
ing current is not gauge invariant (as it would be in the C/(l) case i.e. electromagnetism) and thus 
a function that is not measurable [ 122] - An important property in the SO case is that the gauge 
field is determined by the physical electric field, and thus the gauge is fixed and the spin current 
becomes properly defined. The third term in Eq. I12[ which has the structure of a gauge symmetry 
breaking term, would also by itself fix the gauge and it has further consequences to be discussed 
below. 

Rewriting the last term in Eq ll2l in terms of ordinary derivatives plus a gauge dependent term 
we obtain 

2m 

The second term is the non-Abelian analog of the London term in superconductivity. The main 
result of this approach is then to recognize that such second term exactly cancels the symmetry 
breaking term in Eqll2l and renders zero matter currents proportional to the electric field (zero 
spin conductivity in arbitrary space dimension) . The scenario is now very different from supercon- 
ductivity: there, the London term is the only one remaining after symmetry breaking, while for the 
non-Abelian case, the London contribution gets cancelled. As discussed in references [[TOJ [2T] , equi- 
librium currents remain in relation to the leftover radiative contribution, cubic in the non-Abelian 
potential plus a field independent matter contribution. 

4. Abelian analogy 

In order to discuss an analogy, we will consider the simpler case of the U(l) Abelian gauge 
theory. Then, the Lagrangian density reads as [ [55] 

C = ihip*d t ^ - —(-iKDipYf-iKDrp) - i/>* (e<f>U + -^—B 2 - \e Q E 2 (16) 
2m 2/io 2 

and minimization of the action with respect to the gauge field leads to 

^WVf - f V^) - -If ^ - - V x B = (17) 
2m m fio 

from which the charge current density follows 

e • • e 2 - 

j = — (ip*(-ihVip) + H?iVV0*V0 — M*i>- (is) 

2m m 

The first term is usually referred to as the paramagnetic term, while the second term is responsible 
for diamagnetic properties of matter. If one would have changed the Lagrangian in eq. [11)1 adding 
a GSB term as 

C—>£ + £-Pri>, (19) 
2m 

the current would have been changed to its paramagnetic contribution only, i.e. without any 
dependence on the gauge field (no second term in the r.h.s. of equation [TB]) . This is what occurs, 
as a first approximation, in a paramagnetic metal in a weak magnetic field. The cancellation of 
this diamagnetic term corresponds to the similar scenario of the vanishing of any spin Hall effect 
in the Pauli SO case. 



5. Some implications 

5.1. Vanishing of the spin Hall conductivity for Rashba materials in arbitrary dimension 

As we have mentioned, the presence of the gauge symmetry breaking term in the Lagrangian 
density exactly compensates the "diamagnetic" (also called diacolor) term in the spin current 
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density, and therefore there cannot be any spin current proportional to the electric field [ [5]. 
This means that the spin Hall conductivity identically vanishes, and this result is true in arbitrary 
dimensions. In the particular case of a two-dimensional system with Rashba SO interaction, Rashba 
has shown, using sum rule arguments, that there is no spin Hall conductivity [ 124] . Our conclusion is 
more general in the sense that in any dimension we obtain a spin Hall conductivity which vanishes 
due to an exact cancellation between two terms. In two dimensions, the situation is special in the 
sense that both terms do in fact vanish already! 



5.2. Voltage quantization 

We now consider a two-dimensional electron gas (2DEG) is a crystal, and we analyze the 
consequence of periodicity in the real space. If we consider the transport of the spinor [[25] along 
the primitive cell of the crystal, we have to form quantities like 

exp (iTT 1 <j)(p + gW a T a )dr^j , (20) 

in the absence of magnetic field. If we restrict ourselves to transport along vectors a and b (or 
along b, then along a) in a uniform external electric field, the commutator [Tg, Tg] should be con- 
sidered, where — exp (i'hr 1 (p + gW a T a )dj . This commutator can be shown to be proportional 

to sin 2 ^ '2 \a x E\ sin £ ^ \ b x E\, thus [T^, Xg] vanishes when E is in the plane of the 2DEG [[S]. 
On the other hand, when E is perpendicular to the plane, a quantization condition appears for the 
voltage along at least an arm of the elementary cell, say the arm a 

E z a = pirmc 2 /e (21) 

with p an integer. The quantity 2ir / (2mc 2 / e) plays the role of a quantum of voltage similar to the 
flux quantum in the Aharonov-Bohm (AB) effect [[5]. The analogy is nevertheless not complete, 
and there are important differences. While the gauge invariant phase in the AB effect is given 
by § Adr along a closed path, in the present situation the corresponding quantity reduces to the 
integral along an open path, J \E x dr\. This is of course due to the fact that the non-Abelian 
gauge field is given by the electric field itself, which is gauge invariant already (and there is no 
need to close the path to render the phase gauge invariant) . The present case is thus more similar 
to the Aharonov-Casher effect [126], 

Let us mention that in a semi-conductor, using the effective mass of the electron instead of its 
bare mass, the quantum of voltage would be considerably reduced. 



5.3. Spin interferometry 

In this section we consider an electronic Mach-Zchndcr interferometer where electron beams 
can interfere and are then collected in two distinct detectors [ [37] . A magnetic field perpendicular 
to the plane of the interferometer creates a gauge vector A in the whole space. The magnetic field 
could be limited to a narrow area and does not "touch" the arms of the interferometer. In the 
illustrative case below we consider that the electron mirrors and beam splitters, realized by gate 
potentials, are diagonal in spin space, i.e. they do not mix spin components. Consideration of the 
changes in the electron propagation direction on the spin components [ 28, 31] does not change the 
qualitative scenario drawn below. 

Let us first consider the Aharonov-Bohm situation where the spin of the electrons is neglected, 
and we only discuss the case of one detector. The electrons can follow one of the two paths called 
I and II in Fig. [1] The two electron paths are supposed to be essentially one-dimensional. 

The wave function, transported along each path to detector D a reads as 

i/ia — r 2 exp (ih^ 1 J (p + eA)df^J exp (ih^ 1 j ' (p + eA)df^J t\ 

+t 2 exp(ih~ 1 J (p + eA)dr^J exp (ih^ 1 J (p + eA^dr^r^o- (22) 
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Figure 1. Mach-Zehnder electron interferometer: an electron beam is first split into two distinct 
beams through a Beam Splitter (BS), each of the beams then being reflected on a Mirror (M) 
before being collected at another BS. Two detectors (D) can measure the beams intensity in two 
perpendicular directions. An external magnetic field is applied in the perpendicular direction. 



In order to calculate the probability amplitude in the detector, we have now to evaluate quan- 
tities like 

(23) 



cxp^Ti 1 J Hydi/j exp(^ih 1 J H x dx 



with Hi = pi + eAi (take care to the non standard definition here, in order to shorten the expression 
of the comparator operator [EH])- Using the Baker-Campbell-Hausdorff (BCH) formula and the 
fact that the commutator [7T:c,7ry] is a c— number, this is equivalent to 



exp(ih- 1 (J n y dy + Jn x dx) -^(ih- 1 ) 2 dy J dx[Jly,U x } 



(24) 



The interferences are due to the last term involving [IIy,II x ] = [p y + eA y ,p x + eA x ] — iheB, so 
the last integral contributes to the phase by an amount i7r$s/$o with <&n = h/\e\ the quantum of 
magnetic flux (associated to the charge e). Eventually, one has 

1> a = (r 2 tie^ B ^° +t 2 r ie -^/* )e M Vo, (25) 

where the phase a comes from the non-interfering contribution in Eq. 1241 This result is the stan- 
dard Aharonov-Bohm effect which states that the interference pattern is determined by the total 
magnetic flux enclosed between the two arms of the interferometer. 

We now consider a variant of this problem where the spin-orbit interaction is taken into account. 
The arms of the interferometer are supposed to be made of a "Rashba-Dresselhaus" active medium [ 

ESj. 

Instead of transporting a wave function with a phase variation induced by the coupling of the 
electron charge to the Abelian gauge vector, we now consider the transport of a Pauli spinor and 
its precession due to the coupling of the electron spin to the non- Abelian gauge field. In the case 
when both Rashba and Dresselhaus SO interactions are present, the Hamiltonian is given by the 
following expression [ ESI [3U] , 

H = ^- + V + a(n x (Ty - TT y cr x ) + 0{-K v a y - n x a x ), (26) 
Zm 

where if = p — eA. The non- Abelian gauge field now takes the form 

gW a r a = (/3t x - otT y )u x + (ar x - ^T y )u y (27) 
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and the spinor transport along the paths of the interferometer is defined by the operator 

*a = W«*o, (28) 
according to (here we disregard the fact that the spin components are mixed at each reflection [ 

EH) 

U a = r 2 exp (iTT 1 (n + gW a T a )drj exp (iTT 1 J + gW a T a )dr)h 

+t 2 exp(ih- 1 (U + gW a T a )dr^ exp (ih^ 1 J 0. + gW a T a )drj ri . (29) 

Again, the evaluation of the previous expressions is made delicate due to non commutativity of 
the gauge field components (and their commutator is no longer a simple c— number), so the BCH 
formula is now of no use, since it would require the complete nested expression. We rather use 
properties of the Pauli marices to decompose the exponential functions, exp(±i7<r") = COS7I ± 
ia n sin 7. After some algebra, we get 

U a = .4+ [cos 2 A - sin 2 A sin 20]1 2X 2 + «(sin A)lM2x2, 

where we have introduced the traceless matrix 

1M 2X2 = A- (sin A cos 26)cr z - A+ cos A(cos + sin 9)(cr x - cr y ), (30) 

the dimensionless variables 

A = (m*L/h)y/c^+W, (31) 

9 = tern- 1 (P/ a), (32) 

and the coefficients 

A± = t 1 t 2 ±r 1 r 2 e 2i ^ B ^°. (33) 

In the definition of A, the parameter L is the length of a single arm of the interferometer (the 
two arms are chosen of equal lengths). If the parameter A vanishes, the SO interaction simply 
disappears and we are led to the Abelian AB situation. If cos A = 0, the operator U a is diagonal 
in the input spinor basis, which means that a perfect spin filtering is possible in this original basis 
if one of the two eigenvalues Aj_ vanishes (in which case the corresponding component is filtered 
and only the other component survives in detector a). If it is not the case, filtering may still be 
possible, but in a different basis, tilted from the original one [125]. 

The traceless condition simplifies the diagonalization of IM2X2, and the eigenvalues for U a are 
easily found to be 

A j_ = A+ [cos 2 A - sin 2 A sin 20] =F i sin A x 

^A 2 _ sin 2 A cos 2 29 + 2A\ cos 2 A(l + sin 20). (34) 

The filtering condition (e.g. of the + component of the spinor) is guaranteed by the condition 
A+ = 0. It can be shown that this condition can always be fullfiled by a convenient choice of the 
magnetic flux between the arms of the mterferomter. As we have discussed above, an interesting 
feature is that, depending on the value of A (remember that the Rashba amplitude can be tuned 
experimentally), it is possible to achieve this condition in a non tilted basis. When the perfect 
filtering condition is satisfied, the amplitude of the non filtered component is obtained via the other 
eigenvalue, A^^/g - . In the following figures we propose illustrations in the most generic case, i.e. the 
tilted basis, since this is the situation that occurs in general. In contrast, the condition for perfect 
filtering in the original basis requires a particular relation between the Rashba and Dresselhaus 
amplitudes through the relation cos A = 0. We thus show contour plots of the magnetic flux which 
allows for perfect filtering and for the intensity of the non filtered component in the detector in 
the tilted basis, in the plane a, (3. 
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Weak 

magnetic 

field 



-6 -4 -2 2 4 6 

P 

Figure 2. Perfect filtering by interference for the tilted axis. The plot shows a contour plot of 
sin -k^b/^o in the plane a, f3 (in units of ft/ (m*L)). The darker regions indicate larger values for 
the magnetic flux needed to yield the condition of perfect filtering, from an unpolarized input. 
Highlighted circle and diagonal lines depict the zero flux solutions that yield perfect filtering. 

6. Conclusion 

We have presented a non Abelian gauge theory suited to deal with non relativistic quantum 
mechanics in the presence of various types of spin-orbit interactions. This formulation has the 
advantage of a correct definition of the spin current density and is adapted to treat problems like 
topological quantization and spin filtering. This is a very elegant way to study spintronics. 
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